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Abstract 



We consider the renormalisation of the Fayet-Iliopoulos D-term in a 
softly-broken abelian supersymmetric theory, and calculate the associated 
/5-function through three loops. We show that there exists (at least through 
three loops) a renormalisation group invariant trajectory for the coefficient 
of the D-term, corresponding to the conformal anomaly solution for the soft 
masses and couplings. 



I. INTRODUCTION 



In abelian gauge theories with N — 1 supersymmetry there exists a possible invariant 
that is not allowed in the non-abelian case: the Fayet-Iliopoulos D-term, 

L = £ Jv(x,6,6)d 4 6 = £D(x). (1.1) 

In this paper we discuss the renormalisation of £ in the presence of the standard soft 
supersymmetry-breaking terms 

L SB = (m 2 )&Vj + Q^'V^fc + ^Ffafa + ^MXX + h.c.) (1.2) 

Let us begin by reviewing the position when there is no supersymmetry-breaking, i.e. 
for Lsb = 0. Many years ago, Fischler et al [I]] proved an important result concerning 
the renormalisation of the D-term (see also Ref. p|). Since it is a / d A 8-tjpe term, one 
may expect that the D-term will undergo renormalisation in general. Moreover, by simple 
power-counting it is easy to show that the said renormalisation is in general quadratically 
divergent. Evidently this poses a naturalness problem since (if present) it would introduce 
the cut-off mass scale into the scalar potential. At the one loop level it is easy to show 
that the simple condition Tr^y = (where y is the U\ hypercharge and the trace is taken 
over the chiral supermultiplets) removes the divergence. Remarkably, although one may 
of course easily draw individual diagrams proportional (for example) to Tr3^ 5 , y 7 ■ ■ ■ etc., 
this condition suffices to all orders. 

In the presence of supersymmetry breaking, however, it is clear that £ will suffer 
logarithmic divergences. If calculations are done in the component formalism with D 
eliminated by means of its equation of motion, then these divergences are manifested 
via contributions to the /3-function for m 2 . It is in this manner that the results for 
the soft /3-functions were given in, for example, Ref. ||. Here we prefer to consider the 
renormalisation of £ separately; an advantage of this is that it means that the exact results 
for the soft /3-functions presented in Refs. ||- (see also Refs. || |§) apply without 
change to the abelian case. The result for f3g is as follows: 

& = + & (1-3) 

where fig is determined by V-tadpole (or in components D-tadpole) graphs, and is in- 
dependent of £. In the supersymmetric case, we have fig = 0, whereupon Eq. (|1.3|) is 
equivalent to the statement that the D-term, Eq. ( |1 . 1| ) , is unrenormalised. In the pres- 
ence of Eq. (|Op, however, fig depends on m 2 , h and M (it is easy to see that it cannot 
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depend on b) . The main result of this paper is a complete calculation of (3% through three 
loops; it is interesting that the dependence on h and M arises first at this order. (A 
partial calculation was presented in Ref. [ID|.) 



Although in this paper we restrict ourselves to the abelian case, it is evident that a 
D-term can occur with a direct product gauge group (Gi ® G 2 ■ • •) if there is an abelian 
factor: as is the case for the MSSM. In the MSSM context one may treat £ as a free 



parameter at the weak scale JTTJ] , in which case there is no need to know (3^. However, if 
we know £ at gauge unification, then we need {3% to predict £ at low energies. Now in the 
D-uneliminated case it is possible to express all the /3-functions associated with the soft 
supersymmetry-breaking terms given in Eq. ( |1.2| ) in terms of the gauge /3-function f3 g , 
the chiral supermultiplet anomalous dimension 7 and a certain function X which appears 
only in (3 m 2] moreover in a special renormalisation scheme (the NSVZ scheme), (3 g can 
also be expressed in terms of 7, and X takes a particularly simple form |I2"|| . It is 
clearly of interest to ask whether an analogous exact expression exists for (3^. Moreover, 
there exists an exact solution to the soft RG equations for m 2 , M and h corresponding 
to the case when all the supersymmetry-breaking arises from the conformal anomaly jn| 
and it is also interesting to ask whether this solution can be extended to the non-zero £ 
case. 

The key to the derivation of the exact results for the soft /5-functions is the spurion 
formalism. The obstacle to deriving an analogous result for (3% is the fact that individual 
superspace diagrams are (as already mentioned) quadratically divergent. We do, however, 
present a solution for £ related to the conformal anomaly solution, but which must be 
constructed order by order in perturbation theory. 



II. RENORMALISATION AND NON-PROPAGATING FIELDS 



A. The non-supersymmetric case 



This paper is concerned with the renormalisation of the coefficient of an auxiliary field 
term, and it is perhaps useful to begin with a (we hope) pedagogical discussion of this in 
a non-supersymmetric context. One often sees the statement that the field theory 

L = \{d,<py+ l -F*+ l -hF<f (2.1) 
where <p 2 = Y,a=i a a ; is equivalent to the theory 
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L = \{d^f - (2.2) 

(where 4 = (0 2 ) 2 ), by virtue of the equation of motion for the non-propagating field F, 
which is 

F = - l -h<f (2.3) 

so that 

A' = 3h 2 . (2.4) 

There is a trap for the unwary here, however, in that Eq. ( |2.1| ) is not multiplicatively 



renormalisable, and as a consequence Eq. (|2.4j) is not renormalisation group (RG) invari- 
ant. Let us replace Eq. ( |2.1| ) by 



L = \{d,<P) 2 + V + - A 4_ (2 5) 

We then obtain (eliminating F) 

A' = A + 3h 2 , (2.6) 

and it follows that 

fa =Px + Gh(3 hj (2.7) 

which is easy to verify at one loop by direct calculation: 

16^, = ^±1 X ' 2 (2.8a) 
o 

16vr 2 /? A = ^-^A 2 + 12Xh 2 + 12/i 4 (2.8b) 

lQ^(3 h = ^±i/i 3 + ^±4a (2.8c) 



and it is easy to see that Eq. (|2.7f) indeed holds. The minor subtlety here is that f3\ 
does not vanish when A = 0, so that the naive relation Eq. (|2.4|) is not RG invariant. 
Consequently, if we set A = 0, then Eqs. fl2.8a| ) and ( |2.8cp and are not compatible with 
the (naive) result of taking /i^ of Eq. ( ^4|) . 

One may generalise this example as follows, by introducing a mass for and a linear 
F-term: 
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L = \{d,<PY - \m 2 <? + \f 2 + £F + \hF<? - A 04 

= l(^f-l m 'V-^ (2.9) 

where 

m' 2 = m 2 + h£. (2.10) 

We now have the additional identity 

P m »=P m *+hPt + SP h (2.11) 

which can be verified at one loop using the results: 

2 iV + 2 , 2 

lo7r p m ,2 = — - — A m 

16tt 2 /5 € = ^-h 2 i + Nhm 2 

I6n 2 p m 2 = ^^\m 2 + 2h 2 m 2 (2.12) 

together with the result for p h which is unaffected. 



B. The supersymmetric case: L>-terms 



After this warm-up exercise, let us turn to a softly-broken abelian supersymmetric 
gauge theory. The relevant part of the Lagrangian is as follows: 

L = -D 2 + iD + gD<f>*y<f> - <p*m 2 (f) + ■■■ (2. 13) 

2 

where y i j is the charge matrix of the chiral supermultiplet, and m 2 is a supersymmetry- 
breaking term. After eliminating D this becomes 

L = -<P*m 2 <P- l -g 2 ^y<P) 2 , (2.14) 

where 

m 2 = m 2 + g^y. (2.15) 



RG invariance of this result gives 
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A^m 2 , • • •) = Mm\ ■■■) + My + gfey 

= p m2 (m 2 , ■■■)+ 2(3 9 iy + gy^(m 2 , ■ ■ •), 



(2.16) 



where 
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(2.17) 



with independent of £. For a derivation of Eq. ( |2.17| ), see Ref. [T(J. What about the 
pitfall in the toy model which led us to introduce A? We are saved by supersymmetry: if 
we add a 4 term to Eq. ( [2.13 ), then supersymmetry would be broken, at the dimension 



4 level; contrariwise, if we omit it then it will not be generated. Therefore, Eq. (|2.16|) is 
valid. 



There is an important distinction between fi^ 2 {m 2 



and fi m 2(m 2 



which both 



appear in Eq. ( 2.161) , and determine the mass renormalisation with D eliminated and 
uneliminated respectively. Because /3-functions are determined by 1PI diagrams, fi m 2 
does not contain any D-tadpole contributions; the renormalisation of these is dealt with 
separately by fi^. However, in the D-eliminated formalism, there is no fi%, and there is a 
distinct set of contributions to /3 m 2 involving the four-point vertex created by eliminating 
D. It follows that 



firh,2 (m 2 



(3 m 2 (m 2 



+ ^(m 2 ,---), 



(2.18) 



since diagrams corresponding to one or more insertions of a D-tadpole type contribu- 
tion on the internal line of a diagram do not contribute to the /3-function because the 
corresponding Feynman integral is factorised [ 13 ]. 

Let us now define our notation for the calculation. We take an abelian N — 1 super- 
symmetric gauge theory with superpotential 



and at one loop we have 

16vr 2 /3« = g 3 Q = g 3 Ti [y 2 



167r 2 7 (1)i ,- = P { 



1 -\rikl y 



2g 2 (y 



2\i 



J- 



(2.19) 

(2.20a) 
(2.20b) 



In the spurion formalism the soft-breaking Lagrangian is given by 
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L soit = J d 2 ee 2 {hi iik ^^ h + + ^Mw a w a ) + h.c. 



(2.21) 
(2.22) 



where V is the vector superfield and W a the corresponding field strength. The equivalent 
expression in terms of components is given in Eq. ( |1.2| ). With the explicit all orders result 
for /3 m 2, we prove a remarkably simple result for (3^. The afore-mentioned exact result for 
p m 2 is H 



(/3 m2 ) i j (m 2 , 



where 



r\ r\ «^ 

200* + 2MM*g 2 ^- + Y— + Y*- — + X— 



7' 



0=\Mg 



8 8 

2 17 _ Umn u 



dg 2 



dY 



Imn 



(2.23) 



(2.24) 



Y ijk = (m 2 )\Y ljk + (m 2 ) j iY ilk + (m 2 )^F^ 



and (in the NSVZ scheme) 



167r 2 x NSVZ = _ 2g 3 TT m 2y2 



(2.25) 



(2.26) 



Once again we should emphasise that, whereas in a non-abelian theory Eq. (|2.23|) holds 
in both the D-eliminated and D-uneliminated formalism, in a theory with abelian factors 
it is only true for D uneliminated. 



It is now easy to show that 

/5 m2 (m 2 , • • ■) = Pm^m 2 , ■ ■ •)• 



(2.27) 



This follows simply by substituting for m 2 from Eq. (|2.15|) and then using the facts that 

{y)\Y ljk + (y) j iY ilk + {y) k t Y ijl = o (2.28) 



by gauge invariance, and 



Tr(;y 3 ) = 



(2.29) 



for anomaly cancellation. 
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The result for X, Eq. ( |2.26| ), applies in the NSVZ scheme, which is one of a class 
of schemes related by redefinitions of g and M, the ramifications of which are described 
in Ref. ||. Now X transforms non-trivially under these redefinitions 0, but it can be 
shown using Eqs. ( |2.28| ), ( |2.29| ) that X is unchanged by the replacement m 2 — > m 2 in any 
member of this class of schemes; consequently Eq. (|2.27| ) always applies. We then find 
immediately from Eqs. ( [2.16|) , ( p. 18 ) that: 



k{m\...)=2^ + ^{m\---). (2.30) 

Now on dimensional grounds we may write: 

ft = m 2 A 1 (g, Y, Y*) + hh*A 2 (g, Y, Y*) + MM*A 3 (g, Y, Y*) + (Mh* + M*h)A 4 (g, Y, Y*), 

(2.31) 

where we have suppressed indices for simplicity. (In the conventional DRED 

scheme, will also depend on the e-scalar (mass) 2 , m 2 , and this dependence, as we shall 
see, arises first at three loops. Our three-loop result, therefore, will be in the DRED' 
scheme Hence we have at once that 

Tr(^) = 2% (2.32) 

g 2 

So if we take the .D-tadpole contributions to j3^, then the terms proportional to m 2 will 
reduce to 2/3 g / g if we replace m 2 by g~y. This result is, in fact, clear from a diagrammatic 
point of view, since the aforesaid replacement converts the diagrams into D self-energy 
graphs, and hence indeed gives rise to (3 g . 



III. THE ONE LOOP CALCULATION 



Here we describe the one-loop calculation of this is straightforward, of course. In 
a softly broken theory, the (3% calculation may be carried out in components, or using 
the superfield spurion formalism. Usually, superfield techniques (once mastered) offer a 
substantial reduction in labour compared to component calculations; we will begin to 
suspect from the one loop calculation, however, and confirm in the next section, that this 
is not the case here. 
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Fig. 1: Feynman diagram for the one-loop calculation in components. Dashed lines are 
scalar propagators and the external field is a D. Blobs denote m 2 insertions. 

In components there is a single diagram, shown in Fig. 1, and we have 

= -g(y)'j I d"k (1**, - ^(m 2 )', + • ■ ■) 

f d d k 2 
= -gT±(y) J -p- + gTr(ym 2 ) - + finite terms, (3.1) 

where our integration measure d d k includes the usual (27r)~ d factor. In order to extract 
the ultra-violet divergence from the logarithmically divergent term in Eq. ( |3.1| ) we have 
made the replacement 

d d k r d d k 

(3.2J 



k* J (k 2 + m 2 ) 2 

where mo is an infra-red (IR) regulator mass. Naturally we could have directly evaluated 
the diagram without first expanding in powers of m 2 , but this procedure would be cum- 
bersome at higher loops; it is simpler to treat m 2 as an insertion and introduce regulator 
masses only for those propagators which are IR-dangerous. This technique was described 
in Ref . |1| , and is generally more convenient than the alternative of "threading" a single 



momentum through the diagram |T7| . The pole result for a graph of any number of loops, 
when all sub- divergences are subtracted, is independent of the precise details of how the 
IR divergences are regulated. All this, is of course, well known to higher-loop calculators 
but may, perhaps be of some pedagogical interest. 

We see that to remove the quadratic divergence we must have Tr^y = 0, and that at 
one loop 

k = 7^Tr(ym 2 ). (3.3) 
In the superfield spurion calculation we have two graphs, shown in Fig. 2. 



S 



(a) (b) 

Fig. 2: Feynman diagrams for the one-loop calculation in superspace. Solid lines are chiral 
propagators and the external lines are vector superfields. Blobs denote m? insertions. 

The results are as follows: 

6^(9-9') 



Fig. 2&=-gJ d d k J d i ev(e,e)D 2 [Tr(y) 



k 



2 



+ Tr (ym 2 ) ^D 2 9 2 9 2 D 2 5^(9 - 9')]D 2 \ 0=e , (3.4) 



while 



Fig. 2b = -^Tr (ym 2 ) j d d k j d 4 9 9 2 9 2 V(9, 9)D 2 ^ ^ &) D 2 \ B=e , 

= -gTr (ym 2 ) J d 4 9 9 2 9 2 V(9, 9) J (3.5) 

The first term from Eq. (|3.4|) vanishes via Tr^ = 0, and the second may be reduced using 
the identity Eq. ( |A5| ) to give 

Fig. 2a = sTr (ym 2 ) J ^ J d 4 9V(9, 9)e 26 ^. (3.6) 

If we expand the exponential in Eq. ( |3.6p , the quadratically divergent 9 2 9 2 term cancels 
Fig. 2b, while the remaining term reproduces the component calculation, Eq. (|3.1|). 



IV. THE TWO LOOP CALCULATION 



In this section we discuss the two-loop calculation of (3^ in some detail. Calculations of 
/^-functions for soft-breaking parameters may be carried out in components, or using the 
spur ion formalism. Indeed, as mentioned earlier, in the case of Pm and /3 m 2 the fact 
that the spurion diagrams are only logarithmically divergent means that these quantities 
have simple all-orders expressions in terms of 7 and (3 g . However, as we have emphasised, 
individual diagrams contributing to (3^ are quadratically divergent. This means that if, 
for example, we represent a 7t u vertex in superspace by h l ^ k 9 2 , then we cannot simply 
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factor the 9 2 out, because it can be "hit" by a superspace D-derivative; indeed, as is clear 
from the one-loop calculation, the contribution when the 9 2 is not "hit" will not give a 
logarithmic divergence, and must cancel. The simple relationship between a graph with a 
h^ k and the corresponding one with a supersymmetric Yukawa vertex which holds for the 
soft breaking /3-functions is thereby lost. Nevertheless, the spurion formalism may still 
be used. In this section we shall describe both the spurion approach and the component 
calculation. Normally a superspace calculation would be expected to be more efficient 
than the component version. In this case, however, we shall see that the advantages of 
the spurion calculation are by no means so obvious. The fact that in components the D 
insertion can only be on a scalar line considerably reduces the number of diagrams in this 
case. 




(d) (e) (f) 

Fig. 3: Feynman diagrams for the two-loop calculation in superspace. Solid lines are chiral 
propagators and internal wavy lines are vector propagators. Blobs denote m 2 insertions. 

The two-loop diagrams in the spurion formalism are depicted in Fig. 3. Standard 
superspace manipulations are used to reduce the graphs to basic momentum integrals, 
together with a single remaining / d 4 9; by power counting, the logarithmically divergent 
contributions come from terms with no 9s and 9s remaining in the integrand. Some 
useful identities are collected in Appendix A. Note that we have omitted graphs with a 
mass insertion on the leftmost vertex, where the external V is attached; these graphs, like 
Fig. 2b, do not contribute to the logarithmic divergence, and are cancelled by the quadratic 
divergences (terms with an integrand involving 9 2 9 2 ) from the graphs shown. We have 
also omitted a graph like Fig. 3d, but with the mass insertion on the rightmost vertex, 
because it also gives rise to a quadratic divergence only. The divergent contributions to 
£ B from each graph are listed in Table 1. 
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Fig. 3 


-JS 1 





8JS 2 


-4J5 2 


4JS 2 


-4JS 2 


Fig. 5 


-JS X 


2JS 2 


AJS 2 


-2JS 2 







Table 1: Results for two-loop Feynman diagrams 

Here J denotes the standard two loop momentum integral shown in Fig. 4, and also 



Sx = Tr ym 2 Y 2 



S 2 = g'Tr y 3 m 2 



(4.1) 




Fig. 4 : Momentum integral for the two-loop calculation. The dot denotes a double propa- 
gator. 



The calculation of J proceeds as follows (note that here and in all subsequent integrals 
we subtract all subdivergences): 

d d kd d q 2 



J 



q 2 (k — q) 2 (k 2 + ml) 2 167r 2 e 
2 



drk 



(k 2 + ml) 



2\ 2 



( 4vr )d e 2 

1 

2Y2 



16tt 



1 e \ 4 

H 7e - d 

2 / (47r)l +2 e 2 
2 1 

-- + - 

e z e 



2 



(4.2) 



We have ignored contributions of the form J 2 , where 

d d k 



(k 2 + ml 



2\1 ' 



(4.3) 



such as that from Fig. 3b, because I 2 has no simple pole after sub-divergence subtraction; 
I 2 is the simplest possible example of a factorised Feynman integral, which quite generally 
give no simple pole jyj . Subsequently we will ignore any graph which reduces to factorised 
form. 



Thus using the simple pole given by 

•^simple ~ 



(16tt 2 )V 



(4.4) 
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and recalling that to get the L-loop contribution to the /^-function we need to multiply 
the Feynman diagram simple pole result by L, we find that at two loops we have 



16tt 2 



2gTr ym — 4gTr ym 7 



+ 



(4.5) 



so we see that in fact only A% is non-zero through this order. 



"WV 1 
\ 



'vw' 
\ 



\ 



(b) 



'Wv' 




5: Feynman diagrams for the two-loop calculation in components. Dashed lines 
are scalar propagators, solid lines are chiral fermion propagators, double solid lines are 
gaugino propagators and internal wavy lines are vector propagators. Blobs denote m? 
insertions. 



The calculation may equally well be performed in the component formalism. The 
relevant diagrams are shown in Fig. 5. As we mentioned earlier, there are relatively few; 
indeed, fewer than in the spurion case. Their divergent contributions are again listed 
in Table 1, and upon adding we find again the result of Eq. ( |4.5| ). It is apparent from 
Fig. 5 that there is no DRED /DRED' distinction at this order, because the vector boson 
couples only via the vertex, which projects out the e scalar. A further consistency 

check is provided by Eq. (|2.32j ); since 



16%' P g = g 6 Tr Y - 2g ir Ti ^ 7 



+ 



(4.6) 



we see that Eq. ( [4.5| ) is indeed consistent with Eq. ( |2.32| ). Finally, it is easy to verify that 
our result reproduces the relevant terms from the calculation of f3 m 2 (with D-eliminated) 
presented in Refs. , |L8] . (The other two- loop calculation of the soft /^-functions was 
performed with D uneliminated.) 
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V. THE THREE LOOP RESULTS 



We have calculated ^ 3 ) DRED j n f u n. As we found in the previous section, the cal- 
culation in terms of component fields is generally more straightforward than that using 
the spurion formalism. In the case of terms proportional to w?yY A we have performed 
both spurion and component calculations, while for m 2 y 3 ^ 2 -type terms we have used 
the spurion formalism, which could be streamlined by systematic use of the identities in 
Appendix A. Both these calculations were sensitive to the check provided by Eq. ( J2.32 ). 
The rest of the calculation was done using components. Although the number of diagrams 
is large, the amount of algebra involved in each diagram is not great. 

In both component and superfield formalisms, every graph can be reduced to a sum of 
terms consisting of a product of a group theory factor and one of a set of logarithmically 
divergent three loop graphs, which are shown in Fig. 6. 




D E F 



Fig. 6: Momentum integrals for the three-loop calculation. One/two dots 
denote a double/triple propagator. The two arrows in D denote contracted 
momenta. 



These graphs may be evaluated by the introduction of infra-red regulator masses as 
described for J in the previous section. The results for the simple pole contributions 
(after subtraction of sub divergences) are as follows: 

a 4 1 d 2 1 r 2 i 

^simple 3(167r 2)3 e ' Simple" 3(^2)8^ ° sim P le 3 (167r 2 ) V 

2 1 1 3 1 

AimplC = ~3(167T 2 )V Simple = ^(3)^^^, F simple = - - • (5-1) 



We relegate details of the calculation to Appendix B. The final result may be written 
as follows: 
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167T 



(3)DRED' 



-6(16 



7r 2 ) 2 Tr 



j; m 2 7 (2) _ 4Tr _ r Tr [HH*y] 



+ 2Tr 



P 2 m 2 y -2VC(3)Tr Wy 



+ 12C(3)p 2 Tr M*Hy 3 + c.c. - 144C(3)# 4 MM*Tr y 



(5.2) 



where 



-F 2 m 2 + -m 2 F 2 + K 2 ) l j + 2F w F jpr (m 2 ) r g - 8g 2 MM*(y 2 ) 



>2\i 



^■ = ^F jM + 4/M(y 2 



and 



(167r 2 )V 2 ^ 



lS^ mpi - V(y 2 ) P ^n P n P + 2^ 4 Tr y 2 (y 



,2\i 



(5.3) 
(5.4) 

(5.5) 



with (Y 2 Yj = Y ikl Y jkl , (h 2 Y 3 = a a jkl 



h lkl hjki- We can now check the m 2 terms in this result, 
using Eq. ( |2.32j ). Replacing m 2 by gy, we obtain 



^(16vr 2 ) 3 Tr(yAi° j ) = 6X 1 + 12X 3 + 2X 4 - 12# e Tr y 2 Tr y 



| (3) 



(5.6) 



where 



x 1 




x 3 


= /Tr 


py A 


x. 


= 9 2 Tr 


p 2 y 



knpi 



(5.7) 

in precise agreement with the result for (3^\ given in 119 1, which for an abelian theory is 
(167r 2 ) 3 /?f) DRED = g {3X X + 6X 3 + X 4 - 6/Tr [y 2 ] Tr [y 4 ] } . (5.8) 



(3)DRED 



(Of course for (3 g there is no distinction between DRED and DRED'). Note that (3\ 
would only differ from ^ 3 ) DRED D y inclusion of terms of the form g 5 m 2 Tr [y 5 ] and 
g 3 rh 2 Tr [Py 3 ], arising from e-scalar mass insertions. We have not calculated these explic- 
itly because it is clear they can be removed by a redefinition of m 2 , as follows: 



Sm 2 = - 2 J-rh 2 y 2 + ai U- m 2 y A + a 2 



16tt 2 



167T 5 



m py , 



(16tt 2 ) 



(5.9) 
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where the first term was derived in |15|] . It would be interesting to verify that the appro- 
priate redefinition also renders the three-loop contribution to f3 m 2 independent of m 2 . 

Finally, let us compare our result with the form of ^ 3 ) DRED that we obtained in 



Ref. [|nj (note that we did not there distinguish DRED from DRED' ). We see that our 
result Eq. ( |5.2|) indeed confirms the conjectured form given in Eq. (4.10) of Ref. jnj, and 
that the two then undetermined constants are given by V\ = 24^(3) and ^2 = 0. 



VI. THE CONFORMAL ANOMALY TRAJECTORY 



The following set of equations provide an exact solution to the renormalisation group 
equations for M, h and m 2 : 
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M 

h ijk = _m /3$ 



M ^, 
9 



im 



2\i 



(6.1a) 
(6.1b) 
(6.1c) 



Moreover, these solutions indeed hold if the only source of supersymmetry breaking is the 
conformal anomaly, when M is in fact the gravitino mass. 

This set of soft breakings has generated considerable interest; but there are clear 
difficulties for the MSSM, since it is easy to see that sleptons are predicted to have 
negative (mass) 2 . Most studies of this scenario have resolved this dilemma by adding a 
constant TOq, presuming another source of supersymmetry breaking. A non-zero £ alone 
is not an alternative, unfortunately, as is easily seen from Eq. (|2.15|) ; the two selectrons, 
for example, have oppositely-signed hypercharge so one of them at least remains with 
negative (mass) 2 . This stumbling block may be overcome by introducing an extra U% [20] 



2jJ ; for alternative treatments see Refs. |fL3|] , ||22|| . 

It is immediately obvious that, given Eq. (O), there is a RG invariant solution for £ 
through two loops (for (3%) given by: 



16tt 2 £ = #|M | 2 Tr 



y(7-7 2 ) 



(6.2) 



since differentiating with respect to /x and using Eq. flb.lcp leads at once to Eqs. fl2.17| ), 



(|4.5| ). Interestingly, however, this result for £ vanishes at leading and next-to- leading 
order, since one easily demonstrates that 



Tr 



Ml 







(6.3) 
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and 



Tr 



(2)1 



Tr 



yd 



(1)^2 



(6.4) 



It is interesting to ask whether the trajectory can be extended beyond two loops, and 
whether it in fact continues to vanish order by order. We have shown that there is indeed 
a generalisation of Eq. (|6.2|) to at least three loops (for and that at this order the 
result for £ is non-zero. 



Our result is as follows: 

£DRED' 



g\M \ 



16tt 2 )- 4 {-3/i - 12C(3) (h - 2# 6 Tr [y 2 ] Tr [y 5 ])} 



(6.5) 



where 



h = Tr [yp 3 ] - -(yy j Y jH Y mn p m k p n l + 2<? 2 Tr [y 3 p 2 



2^ 4 Tr y 2 Tr y 3 P 



l2 = 9 2 (y 3 y j Y jkl Y lkm P m l + g 2 Ti y A P 2 + 2<? 4 Tr y b P 



(6.6) 



It is easy to verify that the result of taking fj,-^- of Eq. ( |6.5| ) is identical to that obtained 
by substituting Eqs. (|6.1j) in Eqs. (|4.5|), (|5.2|) . This is a non-trivial result in that the 
number of candidate terms for inclusion in Eq. (|6.5|) is considerably less than the number 
of distinct terms which arise when Eq. ( |4.5| ), ( |5.2| ) are placed on the RG trajectory. We 
therefore conjecture that the trajectory extends to all orders. 

~ (3) 

It is natural to ask what the result for (3^ is in the NSVZ scheme, which is obtained 
(at the relevant order) by the redefinitions @] 



(16tt 2 ) 2 ^ = -^ 3 Tr [Py 2 



(16tt 2 ) 2 5M = -Mg 2 {Tr [py 2 ] - 2g 2 Ti [{y 2 ) 2 }} + ~g 2 h ikl Y jkl (y 2 y V 



(6.7) 



It is straightforward to show that in order to obtain the results Eqs. (|2.17 ) and ( 2.32|) in 
the NSVZ scheme, we must also redefine £ as follows: 



(16tt 2 ) 2 <5£ = —g 2 Tr [ p y 2 ] f - £ Tr [m 2 Py 



(6.8) 



The effect of this is to replace Eq. ( |6.5| ) by 

£NSVZ 



g\M \ 



(16tt 2 )- 4 {-4/x - 12C(3) (l 2 - 2g 6 Tr [y 2 ] Tr [y 5 ])} 



(6.9) 



16 



and Eq. (O) by 



A(3)NSVZ 
(l 67r 2 )3 ^ 

9 



-4(16vr 2 ) 2 Tr ym 2 7 (2) - - (2Tr [WPy] + Tr [HH*y]) 



- 24/C(3)Tr ly^ 3 + 12C(3)/Tr M*Hy 3 + c.c 



144C(3)# 4 MM*Tr y 



(6.10) 



It is disappointing that this expression does not immediately suggest an all orders result. 
At this point it is worth recalling that, while to connect the DRED' and NSVZ schemes via 
Eq. ( |6.7|) we redefined g and M, there exists also a redefinition of Y (involving C(3)) which 
has the pleasant property of extending to three loops the existence of finite N = 1 theories 
Unfortunately this redefinition disturbs Eq. ( |2.23 ), which leads one to imagine 



that there might be a combined redefinition of m 2 , Y that both preserves Eq. (|2.23|) and 



> (3) 

simplifies (3^ . We have not yet succeeded in constructing such a transformation. 



VII. FINAL REMARKS 



We have presented a detailed, and we hope a reasonably self-contained description 
of the calculation of through three loops. It is intriguing that in the abelian case we 
are unable to express the renormalisation of the theory completely in terms of j3 g and 
7, which, in the non-abelian case, suffice to describe the renormalisation of both the 
unbroken theory, and also the theory with the standard soft terms. Although there exists 
perturbatively a solution related to the AMSB solution for the soft parameters, once again 
we are unable at the moment to extend this solution to all orders. 

The next step is obviously an extension of our calculation to the case of a product 
gauge group including both abelian and non-abelian factors, such as the MSSM; this is 
not a trivial deduction from the results we have presented. Although it is clear that if 
£ is assumed to be small at gauge unification, then it does not have much effect at low 
energies, it should be remembered that this is an assumption, and that the MSSM has 
one more parameter than is commonly supposed. 
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APPENDIX A: D- ALGEBRA IDENTITIES 



In this Appendix we list some identities that we found useful in superspace calculations 
of contributions to (3^. (An early reference for superspace calculations incorporating soft 
breaking is Ref. ||24j| .) The soft terms given in Eq. ( p.22|) are treated as insertions in 
the superfield diagrams and standard superspace manipulations may then be used to 
reduce supergraphs to ordinary momentum space integrals. It is convenient to denote the 
momentum space version of the superspace covariant derivatives by 



{Dpi) a 

(Dpi) A 



d 



86$ 



yaoPl 



d 

de? 



yaa u l 



(Al) 



where 



aa : M' 



(A2) 



We then have the fundamental supersymmetry algebra 

{(Dpl) a , (D p l)a] = \i>oca- 



(A3) 



We find 



aV^-e^ 1 , D 2 Jl = -e 



(A4) 



and moreover 



-{e 1 p 1 +d2p2) 



(A5) 



where 



8 12 = 5^ (6i - 6 2 ) . 



(A6) 



We also have 



^D 2 ql D 2 ql e 2d ^D 2 ql D 2 ql 5i 2 



5i2 

812 



[p - q f - 2 P 2 e 1 0i + 2q 2 e 1 p 1 + p 2 q 2 e\el 

(p + q) 2 + 2p 2 6i0 1 + 2q 2 9 1 1 + p 2 q 2 9 2 1 9 2 1 



(A7a) 
(A7b) 
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Finally, 



5 12 D 2 rl D 2 rl [(p " qY - 2p 2 Oi0i + 2q%p i +p 2 q 2 6 2 1 6 2 1 ]D 2 1 D 2 . 1 5 12 

— 812 (p — q) 2 r 2 + 2q 2 p.r — 2p 2 q.r + p 2 q 2 + terms in 9\,9~i, (A8a) 
St^D^ip - q) 2 - 2p 2 6 l l + 2q 2 6 l f6 l + p 2 q 2 9 2 9 2 }D 2 1 D 2 1 5 1 2 

— 8x2 (p — q) 2 r 2 — 2q 2 p.r + 2p 2 q.r + p 2 q 2 + terms in (A8b) 



Note that the right-hand sides of Eqs. ( |A7a| ), ( |A7b| ) are related by q <-> — q, and similarly 
those of Eqs. (|A8a|) , (|A8b|) are related by r <-> — r. 



APPENDIX B: THREE LOOP DETAILS 



In this Appendix we give a complete graph-by-graph description of the three-loop 
calculation. We start by giving a list of the distinct tensor structures involved: 



T\ = (Y 2 y^ kl Y lkm (m 2 y) m i, T 2 = (Y ) l jY 3 Yi 



T 4 = g 2 Tr \Y 2 m 2 y 3 
Ts = g 2 Y m Y imv Vn> 



T 5 = g 2 Y ikl Y imn (m 2 y) m k (y 



[m 



2\n 



2\m -\ in 



y 



Tq = Tr 



Y 2 Y 2 m 2 y 



1 ...2 \m/-yZyn^ r ,, l.-p,. ,,,*- v"> 



T 7 = £ 4 Tr m 2 ^ 5 , T 8 = g A Ti y 



Tr 



T 9 = Y m Y imn hjj t ih pmn y p T 10 = Tr[Y 2 h 2 y], T n = g 2 Tr[h 2 y% 

i 

3i 



T 12 = ^M/^F^y 3 )^ T 13 = ^ 2 MM*Tr[F 2 ^ 3 ], T 14 = g 4 MM*Tr[y 5 } 



(Bl) 



We now give a list of diagrams contributing to these tensor structures, with the exception 
of T7, for which there are a very large number of separate diagrams; note that when 
m 2 is replaced by gy, T 7 and only T 7 produces Tr [3^ 6 ], so that we can in fact infer the 
coefficient of T7 in our final result via Eq. ( |2.32| ). We did, however, perform the explicit 
T 7 calculation, and indeed obtained the expected result. 
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(d) 

Fig. 7: Feynman diagrams in superspace for the three-loop contribution of the form 
m 2 Y A y , i.e. T1...3. Blobs denote m 2 insertions. 




Fig. 8: Feynman diagrams in components for the three-loop contribution of the form 
m 2 Y 4 y, i.e. T1...3. 
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A > N 
/ \ / \ 

\ / ; 

\ \ / ' 

\ \ / / 

(a) 



(b) 



(C) 



\ ' 

\ / 



(d) 

Fig. 11: Feynman diagrams in components for the three-loop contribution to T 9 ,T 10 . 



-vwl 



(a) 



(b) 



(c) 



(d) 



ax 



(e) 



\ \ 



(f) 



-wv) K/vwi 



•VW| 



(g) (h) (i) 

Fig. 12: Feynman diagrams in components for the three-loop contribution to Tn. 



(a) 



(b) 



Fzg. 15 1 .- Feynman diagrams in components for the three-loop contribution to T± 2 . Blobs 
denote gaugino mass insertions. 
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(a) 



(b) 



(c) 



'\/w{ 



(d) 

1^: Feynman diagrams in components for the three-loop contribution to T 13 . 



(a) 



'wv! 



(b) 




(c) 




(d) 



(e) 



(f) 




(g) 



(h) 



Fig. 15: Feynman diagrams in components for the three-loop contribution to T i4 . 

We begin with a comparison between superspace and component formalisms. The 
results from Fig. 7 (the superspace calculation) are 

Fig. 7a = - 1 -{B + 2L>)(T 1 + T 2 ), Fig. 7b = A7\, Fig. 7c = l -CT 2) Fig. 7d = ^BT 3 ; 

(B2) 

while from Fig. 8 (the component calculation) we find: 

Fig. 8a = -\{B + C + 2D)(T! + T 2 ), Fig. 8b = (A - F)T U 

Fig. 8c = FT U Fig. 8d = ]p{T x + 2T 2 ), Fig. 8e = ^BT 3 . (B3) 
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(Here and elsewhere, we combine diagrams which clearly give identical results, up to 
symmetry factors. Consider for example, Fig. 8d. Because the scalar fields are complex, 
this diagram represents two distinct (by the usual rules) Feynman diagrams.) The totals 
of Figs. 7, 8 are manifestly identical, which is a good check on our spurion rules. 

The results from Figs. 9, 10, 12, 14, 15 are given in Tables 2-6 respectively: 



a 


-2E 




b 


-AA 


\Ta + \T 5 + T 6 


c 


-AA 


T 5 


d 


-2C 


-T 5 - T 6 


e 


-2(2,4 + B + 2D) 


1 n~> 
-2^4 


f 


-AA 


-i t 4 + |t 5 


g 


AA 


1 T 1 

— 2 4 


h 


AA 


1 ' i ' 1 / / 1 ' / 1 
4 J 4 — 2 J 5 ~ J 6 


i 


-4(44 — C — E) 


1 Tl 

-2^4 


j 


-A{2A -2D - E) 


1 ■' / 1 1 ' / 1 ' / 1 

4J4 — 2 Js _ J 6 


k 


-8(A + D) 


1 ■' / 1 1 ' / 1 ' / 1 

4J4 — 2 J 5 — J 6 


1 


-AA 


1 Tl 

-2^4 


m 


-AA 


1 T 1 

-2^4 


n 


2B 


T 4 


o 


2A 


T 4 


P 


2B 


T 4 


q 


AA 


1 1 T" 1 T 1 
— 2 i5 _ ie 


r 


AA 


1 rri 


s 


—A(A + B) 


T 4 


t 


-2(2A + B + 2D) 


T 6 


u 


2A 


T 4 


V 


-2A 


T 4 


w 


-2B 


T 4 


X 


—A(A + B) 


T 4 


y 


-2B 


T 4 


z 


2B 


T 4 


aa 


2B 


T 4 


bb 


-2B 


T 4 



Table 2: Results for Fig. 9. 
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a + b 


c 


d +e 


f 


g 


Fig. 10 


(C-4A- 2B) 


A(A - F) 


(25 + 2F - AA) 


2F 


C 



Table 3: Results for Fig. 10 (all multiplied by T$) 
The results from Fig. 11 are given by 

Fig. 11a = ~AT 9i Fig. lib = ^BT 9 , Fig. 11c = ~\aT w , Fig. lid = ^BT W . 

(B4) 





a 


b 


c 


d 


e 


f 


g 


h 


i 


Fig. 12 


-A-B 


.4 


-2B 


-A 


2B 


1* 


—A 


\{AA-B - 2E) 


A-B-E 



Table 4- Results for Fig. 12 (all multiplied by T lly ) 
The results from Fig. 13 are: 

Fig. 13a = (E — 2A)(T 12 + (T 12 )*), Fig. 13b = B(T 12 + (T 12 )*). (B5) 





a 


b 


c 


d 


Fig. 14 


-C 


-(A + B + 2D) 


2C 


2{A-B) 



Table 5: Results for Fig. 14 (all multiplied by T 13 ) 





a 


b 


c 


d 


e 


f 


Fig. 15 


%(E + \C-A + \B) 


4(5 - C) 


-2C 


8(5 - A) 


8(i£ + |5- A) 


-2C 




g 


h 










Fig. 15 


4(A + B + 2D) 


8(A-5) 











Table 6: Results for Fig. 15 (all multiplied by Tn) 

The final total is obtained by combining the tables, substituting the simple pole results 
for A,B---F from Eq. ( |5.1|) , and multiplying by 3 (for 3 loops): 
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A(3)DRED' 

(16tt 2 ) 3 -L = TTi + 4T 2 - |T 3 + [10 - 24C(3)] T 4 - 12T 5 

+ 16 [1 - 3C(3)] T 6 - 16T 7 - 12T 8 - §T 9 

- 2T 10 + 8 [1 - 3C(3)] T n - 2 [5 - 6C(3)] (T 12 + T* 2 ) 

+ 16T 13 - 8 [13 - 18C(3)] T 14 , (B6) 

which can easily be recast into the form given in Eq. (]5.2|) . (As indicated earlier, we have 
suppressed details of the T 7 computation). 
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